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Abstract. I investigate the possibility of condensation in ideal Fermi systems of 
general single particle density of states. For this I calculate the probability wpf Q 
of having exactly Nq particles in the condensate and analyze its maxima. The 
existence of such maxima at macroscopic values of No indicates a condensate. 
An interesting situation occurs for example in ID systems, where uijv may have 
two maxima. One is at No = and another one may exist at finite No (for 
temperatures bellow a certain condensation temperature). This suggests the 
existence of a first order phase transition. The calculation of wn allows for 
the exploration of ensemble equivalence of Fermi systems from a new perspective. 
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1. Introduction 

Huge progress have been made in the last decade in cooling trapped Fermi or Bose 
gases down to the Fermi degeneracy temperature or to the onset of the Bose-Einstein 
condensation, respectively (see for example Refs. ^ HI El EH an d citations therein). 
Tuning experimental parameters, like the trap frequency or the magnetic fields applied, 
one can change the strength of the interaction between the particles and the density 
of the gas. This freedom, which allows for more accurate comparison between theory 
and experiment, led to a burst of scientific research in the field. The quantum 
effects in these gases (fermionic degeneracy and Bose-Einstein condensation) are 
usually observed at temperatures around or below 1 /iK. Such temperatures can only 
be achieved by evaporative cooling [5], but the time scale for the cooling process 
depends strongly on the elastic collision rate of the particles, which determines the 
thermalization rate [H]- In the case of fermions, at temperatures below the Fermi 
temperature, the scattering probability falls dramatically due to the Pauli blocking 
mechanism, which makes cooling very slow. In contrast to this, in the case of bosons, 
the presence of the Bose-Einstein condensate enhances the cooling rate due to the 
assisted scattering of hot particles out of the trap, and the cooling process can continue 
until almost all the gas is condensed. In analogy to the bosonic case, the existence of 
a condensate in the Fermi system would create the possibility of producing a highly 
degenerate Fermi gas by removing (evaporating) the particles above the condensate. 
Moreover, the Fermi condensate might have the same role of enhancing the cooling 
rate, as the Bose-Einstein condensate. It is well known that a Fermi gas with BCS 
interaction undergoes a phase transition to a condensed, superconducting phase 
In this paper I shall show that a condensate may exist even in ideal Fermi gases. 

At low temperatures, there will be a number No of fermions that occupy the lowest 
single-particle states. If I number the single particle states from to oo, starting from 
lowest energy level and going up-wards, then the first unoccupied state is numbered 
Nq + 1. (The order of the degenerate levels is not important for macroscopic systems.) 
If -/Vq is a macroscopic number, then I will say that it forms a condensate (see figure 



Usually, at low temperatures the particles below the Fermi level are said to be on 
the "Fermi ground state" (see for example Refs. 8 ). The only distinctive physical 
property of the Fermi energy is that at (very) low temperatures, the particle and hole 
populations are symmetric with respect to it. Due to the analogy between the Bose- 
Einstein condensate (see Ref. (Hj) and the Fermi condensate I suggest that the Fermi 
condensate may be better suited for the role of Fermi ground state. 

In Appendix B of Ref. JU| I showed that a condensate forms in an interacting 
system (with constant density of states), with general microscopic exclusion statistics 
properties At the condensation temperature the system undergoes a first 

order phase transition. In I showed how the condensation in the corresponding 
nonintcracting system changes the character of this phase transition. 

2. Condensation 

Let me consider a Fermi system at temperature T and chemical potential fi. Its 
grandcanonical partition function is 




(1) 



in 
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Figure 1. Single particle energy axis. Single-particle states are numbered from 
to oo, starting from the lowest in energy and going up-wards. 



where the sum is taken over all the microstates to, of internal energy U m and particle 
number N m ; as usual, (5 = {k^T)^ 1 . If the system is ideal, then N m = Y^aLq n i,m. an d 
U m = J2i*Lo n i.m,£i, where by e, I denote the energy of the single particle states. I take 
e = and < e i+1 , for any i > 0, as mentioned in the Introduction section. On each 
single particle state there are ni im = or 1 particles. Introducing the expressions for 
U rn and N m into equation I can rewrite Z in the well-known form, 



2 = IIE cxph/3^ - n)m] = J] i 1 + exp[-/3(e, - m)]} • (2) 

i = 0m=0 i=0 

The logarithm of Z is related to the grandcanonical thermodynamic potential, 



f3Q = logZ = lQ g l 1 + e 0ifl - u) 



(3) 



The probability of the microstate to is w m — Z~ l exp [— (3(U m — [iN m )]. 

I assume that for particles in an arbitrary trap and arbitrary number of 
dimensions, the density of states (DOS) has the expression a(e) = Ce s , where C 
and s are constants (I take s > —1). For example if particles of mass to are in a 
d-dimensional box of volume V, then C = V(2Tr)~ d / 2 T- 1 (d/2 + l)d{m/h 2 ) d / 2 and 
s = d/2 — 1. I shall use the same procedure as in Ref. 10 , to calculate the number of 
condensed particles. For this, let me calculate the probability that the lowest iVo + 1 
states are completely occupied, i.e. no<i<Ar = 1 and rii>jv may be either or 1. The 
probability of such a configuration is 



w Nq+ i = Z ■ exp 



• exp 



N 



PKNo + 1) 



i=N + lni=0 



[j {l + exp[-/3( ei -/i)]} 



(4) 
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Using the expression for a to transform the summation Xa=o e * an integral and 
denoting Z ex (N , (3, 0ji) = Yl°^ No+1 + ex P[~P{ € i ~ equation 10} becomes 



WN + 1 



exp 



s+2 



-0C- 



+ (3^C- 



Zck(N ,/3,/3h) _ Z No 
Z ~ Z 



(5) 



where eo is given by the equation Nq = Ce s +1 /(s + 1) (is the level up to which all the 
states are occupied by No + 1 particles), while Z ex (No, f3, 0fj) is the partition function 
of the particles above the level eo- Note that J2n Wn o ^ I an d Zn =o = Z. wn +i 
is the probability that at least No + 1 particles are condensed. 

The quantity of direct relevance here is the probability to have exactly No particle 
condensed. Such a configuration is obtained by removing the particle from the level 
£jv — e o- Doing this in equation J3J or JS}, I obtain the probability 



exp 



-13 [C 



eo + [C 



8+1 



Z ex (N ,f3,Pn) _ Z No 



(6) 



Note that w eQ = cr(eo)wAr ( eo ) is the probability density along the energy axis. 
Nevertheless, in what follows I shall work with wn - If I denote the average of No by 
Nq = ^2n q Now n i then No is macroscopic if (but, in principle, not only if) wn has a 
maximum at No <max , so that in the thermodynamic limit N 0tinax /N > 0, or £o,max > 
[where No,m&x = ^ e omax/( s + !)]• Since Z does not depend on Nq, the maximum of 
wn may be found by solving the equation d\ogZ No /dN a — 0, or d\ogZ Na /de = 0. 
First I calculate log-Zjvb from equation ©: 



log Z. 



,s+2 



c 



s + 2 
de e s log 



eo )+ 0H [C 



r s+l 



1 + e 



(7) 



where the integral represents logZ cx (No, P, /3/i). The derivative of equation (J7J) gives 



d log Z 



N 



de 



= - C[3el +1 + 
log 



CP^Q-Ce^log 



1 + c _/3 ^ 0- ^ 



Ce° 



I _|_ e /3( e o-M) 



JL 

Ce s r , 



or, without assuming anything about a, d\ogZjq /dto = — cr(eo) {log [l + e^ 0- ^] — [er^o^sT] -1 }. 

The interpretation of equation (JHJ) is simple. First note that dlogZN /deo is 
continuous on the interval (0, oo). If eo ~ ^jl, then log [l + e^ 60 ^'] is of the order 1, 
so dlogZ^ol deo < 0, since in the thermodynamic limit [a(eo)kBT}~ 1 <C 1. Moreover, 
for /3(e - n) » 1, 

—dTo Ce ^Y°-^-C4 



and because s > — 1, <91ogZjv /9eo becomes negative for large enough eo, as one would 
expect (the probability to have an infinite number of particles condensed is zero). To 
study in more detail the existence and the number of solutions, let me divide the 
problem into three cases: (1) s — 0, (2) s < 0, and (3) s > 0. 
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Case (1) is the simplest and corresponds, for example, to a Fermi gas in a two- 
dimensional (2D) box. This case may be analyzed in connection with the Bose-Fermi 
thermodynamic equivalence in (2D), which was outlined in Section 2.2 of Ref. |10|. 
The properties of ideal 2D Bose gases have been studied extensively in the past (see 
for example an d citations therein). For s = equation (JSJ) becomes: 

d\ogZ Na 



de Q 



s=0 



If I take the thermodynamic limit simply as {Ck^T)^ 1 = 0, then 91ogZjv /(9eo| s=0 < 
for all eo an d T, so there is no condensate. On the other hand, note that for 
/3(e - £t) < -1, log [l + e^ 60- ^] « e 0^o-iA and, because of the exponential 
dependence on /3(cq — /i), the term in curly brackets in equation 10 may become 
negative even for a macroscopic system. Since log [l + e^ e ° _M )] is monotonically 
increasing with eo, while (Ck^T)^ 1 is a constant, equation JSJ) has maximum one 
solution, and this must be for < eo < A*- I define the condensation temperature 
2c, 2D by the equation {d log Zjv / deo)r=T c 2 d,«o=o = 0- I 11 the approximation /3/i 3> 1 
this condition gives: 

— ^— = io g [(fc B r c , 2D )c]. (io) 

For any T > T C 2D, equation ||5J has no solution, while for T < T" C .2D it has one and 
only one solution for eo > 0, so the system is condensed. To see if equation (fTOf) has 
any relevance for macroscopic systems, let me consider a 2D gas of electrons of Fermi 
energy ep 1 = 11.7 eV (Fermi energy of electrons in Al). If I take the area of the 2D 
gas to be S = 1 m 2 (surely a macroscopic system) and fi = ep 1 , then 

T c A i D « 3263 K. (11) 

At room temperature T r = 300 K, keeping fi — ep 1 , I obtain e^ max « 10.7 eV, only 
about 1 eV smaller than fi. If I take Ao ~ iVo, max , then the condensate fraction 
(Ao/A) is e^U/eAi « 0.91. 

The condensation energy scales roughly with the chemical potential. For another 
example I take heavily doped silicon with the doping concentration of n = 4 x 10 25 m -3 
|21|. For this concentration, the Fermi energy is ep 1 rs 0.04 eV, which gives a 
condensation temperature of about 14 K. In this case, the condensate fraction at, say 
100 mK (roughly the working temperature for the microcoolers or microbolomctcrs), 
is over 99%. 

Two systems of equal number of particles are called thermodynamically equivalent 
if they have the same entropies as function of temperature. All 2D ideal gases of equal 
DOS are thermodynamically equivalent, irrespective of their microscopic exclusion 
statistics [121 HH IS1 CS1 EZl EHI CHI EOl CD] , and this is due to the similarity between 
their excitation spectra |10j . If I put in correspondence the Fermi system under 
investigation with the equivalent Bose system, and I denote by /is the chemical 
potential of the later, then fi^ = fi — ep (where ep is the Fermi energy). The number 
No of fermions in the condensate is equal to the number of bosons on the ground state, 
say -Ab,o E3- Moreover, in the canonical ensemble the condensate fluctuations are 
the same for the two, equivalent, Bose and Fermi gases. 

Case (2) is maybe the most interesting. For s < 0, <9 log Zjv /9eo < in both 
limits, eo — > and eo — * oo, so eventual solutions of the equation <9 log iJ/v /<9eo = 
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come in pairs. To see that there is only one pair of solution note that that 
the derivative Slog [l + e^ 60- ^] /de = /3/[exp(/3/i, — /3e ) + 1] increases with e , 
so log [l + e' 3 ^ 0- '''] is a function concave up-wards, while 1/Ceg is a function 
concave down-wards (—1 < s < 0). Therefore the two curves (log [l + e^ e ° -A1 )] 
and (3/Ccq) may cut each-other either in exactly two points, or in zero points. 
Moreover, if I assume that they do cut and at eo ~ fJ-, or bigger, then k B T ■ 
{<91og [l + e^ 60 ^] /deo} — [exp(/3^i — (3e ) + is of the order one or bigger, while 
-k B T ■ [d(l/Ceg)/de ] = ~sk B T/[Ce s k B T ■ (e /k B T)] < 1, in the thermodynamic 
limit. Therefore, <9{log |"l + e^ 0- ^] — -J^y/deo > for e ~ Mi or bigger. Since 



log [l + e /3(e °"' j) ] - -^g > 0, for eo » /i, I conclude that dlogZ N Jde = has no 
solutions for e G [fi, oo), no matter what is the value, or sign of /i. If 91ogZjv /9e = 
has solutions, then fi > and the solutions are in the interval eo <= (0, /J.)- 
I define the condensation temperature by the equation 

dlog-ZjVo 



max 



= 0. (12) 

T=T C , 3<0 



For T > T c s< o, equation JSJ has no solution, so the system is not condensed. Particles 
on a one-dimensional (ID) interval have er(e) = Ce -1 / 2 , so I shall take ID electron 
systems as examples. If I take again the Fermi energy of Al for electrons on an interval 
of length ! = lm,I obtain T Ci i D » 6242 K. 

The interesting aspect about the s < case is that below the condensation 
temperature Z/v (eo) has two maxima: one at eo = and one at eo = eo,max > 0. By 
comparing the areas below the peaks, one can conclude which maximum is the stable 
solution and which is the metastable one. Similar to the situation outlined in Appendix 
B of Ref. JOj, this leads to a first order phase transition from an uncondensed, to a 
condensed phase, in ID canonical Fermi systems. 



Case (3) includes three-dimensional (3D) ideal systems, which correspond to s = 
1/2. If s > 0, obviously d\ogZ No /de > at e = 0. Since log [l + e^ "^] 
is monotonically increasing and (CfceTep) -1 is monotonically decreasing, while 
<9 log i?7v /9eo is negative at large eo, I conclude that exists one and only one finite 
eo, max, and therefore a finite -/Vo, m ax, at which wn max attains its maximum. For 
example for one cubic meter of Al, at room temperature, eo, m ax ~ 11.2 eV and the 
condensate fraction is about 94 %. 

An important example is taken from the experiments performed by D. Jin and 
collaborators [51 [3] and I shall use parameters from Ref. 0. I describe N = 6 x 10 6 
atoms of 40 K (I take both spin polarizations), confined in a 3D harmonic trap of 
(geometric) mean frequency uj/2-n: = 70 Hz (Huj w 2.9 x 10~ n eV). The DOS is 
o-(e) = e 2 /(?iw) 3 and the Fermi energy is e F = (37V) 1/3 S.u; w 7.6 x 10~ n eV. With this 
parameters I calculate two values of eo !max : £o,max 

(Ti = 1 fj,K) w 2.6 x KT 12 eV and 
eo,max(T 2 - 300 nK) w 1.3 x lO" 11 cV.'Note that both, e , ma x(2i) and e , ma x(T 2 ) are 
smaller than hw, so there is no condensate, although the equation <9 log Zjv /3eo = 
has a solution. Eventually the possibility to obtain a highly degenerate Fermi gas is 
to agglomerate may particles in a trap and then to separate the condensate formed at 
the "bottom" of the trapping potential. 
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3. Conclusions 



In this paper I discussed the formation of a condensate in ideal Fermi systems with the 
density of single particle states of the form a(e) — Ce s (where C and s are constants). 
I did this by calculating the probability wn of having exactly Nq particles in the 
condensate. If s > (like in 3D boxes or traps), wn has a maximum for Nq > 0, which 
indicates the formation of a condensate. The maximum persists at any temperature, 
but eventually N becomes microscopical as T increases. For a cubic meter of Al at 
room temperature, N /N w 0.94 (where N is the total number of particles in the 
system). Unfortunately for the current experiments in harmonic traps, Nq appears to 
be smaller than one, which means that there is no condensate. Maybe a high increase 
of particle number in the trap, without much concern for the temperature, would lead 
to the formation of a condensate, which may be afterwords separated by evaporating 
the uncondensed particles. 

The case s — corresponds to 2D boxes or ID harmonic traps. This is equivalent 
to a 2D Bose gas and the results of canonical calculations may be interchanged down 
to the microscopic scale, with proper redefinition of single particle energies (see Ref. 
|10|1- After the interchange, the particles in the Fermi condensate become the particles 
onto the Bose ground state, which have been extensively studied [T2| . 

Maybe the most interesting case is s < 0. Here wn may have two maxima. One 
maximum is always at Ao = (no condensate), while the other maximum forms at 
finite Nq, for temperatures below a condensation temperature. The existence of the 
two maxima suggest a phase transition of order one in canonical Fermi systems and a 
reconsideration of the ensemble equivalence. In this category are included ID Fermi 
systems. 
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